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ABSTRACT 

We analyse 15 XMM-Newton observations of the Seyfert galaxy NGC 4051 obtained 
over 45 days to determine the ultraviolet (UV) light curve variability characteristics 
and search for correlated UV/X-ray emission. The UV light curve shows variability 
on all time scales, however with lower fractional rms than the 0.2-10 keV X-rays. On 
days- weeks timescales the fractional variability of the UV is F var ~ 8%, and on short 
(~ hours) timescales F vsir ~ 2%. The within-observation excess variance in 4 of the 15 
UV observations was found be much higher than the remaining 11. This was caused by 
large systematic uncertainties in the count rate masking the intrinsic source variance. 
For the four "good" observations we fit an unbroken power-law model to the UV power 
spectra with slope — 2.6±0.5. We compute the UV/X-ray Cross-correlation function for 
the "good" observations and find a correlation of ~ 0.5 at time lag of ~ 3 ks, where the 
UV lags the X-rays. We also compute for the first time the UV/X-ray Cross-spectrum 
in the range 0-28.5 ks, and find a low coherence and an average time lag of ~ 3 ks. 
Combining the 15 XMM-Newton and the Swift observations we compute the DCF 
over ±40 days but are unable to recover a significant correlation. The magnitude and 
direction of the lag estimate from the 4 "good" observations indicates a scenario where 
~ 25 % of the UV variance is caused by thermal reprocessing of the incident X-ray 
emission. 
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1 INTRODUCTION 

The broad-band spectral energy distribution (SED) of Ac- 
tive Galactic Nuclei (AGN) requires multiple emission com- 
ponents, including thermal and non-thermal mechanisms 
l|Shang et al.ll2011T ). The energy output in Seyfert galaxies 
and quasars is dominated by the ultraviolet (UV) emission, 
thought to be mostly thermal emission from the inner parts 
of an accretion disc surrounding the central supermassive 
black hole (SMBH). The location of the UV emitting re- 
gion depends on the details of the accretion flow, and these 
in turn depend on the black hole mass and accretion rate, 
but is typically ~ 10 — 1000> 9 (where r g = GMbh/c 2 is 
the gravitational radius). By contrast the X-ray spectrum 
is usually interpreted as the result of inverse-Compton scat- 
tering of soft thermal photons by an optically thin corona 
of hot electrons in the central few tens of r g from the black 
hole (e.g. Haardt & Maraschi 1991). 

The causal connections between these processes are still 
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unclear but can in principle be investigated by studying the 
time variations in the luminosity across different wavebands. 
Strong UV and X-ray variabil ity is commo n to AGN on a 
wide range of timescales (e.g. ICollinl 2001 ), with the mos t 
rapid variations seen in X-rays fe.g. lMushotzkv et al.|[l993l ) . 
If the emission mechanisms are coupled the UV and X-ray 
variations should be correlated, in which case the direction 
and magnitude of time delays should reveal the causal rela- 
tionship. Two favoured coupling mechanisms are i) Comp- 
ton up-scattering of UV photons — produced in the disc — 
to X-ray energies in the corona (|Haardt fc M araschi 199ll ). 
and ii) thermal reproce ssing in the disc of X- ray photons 
produced in the corona (|Guilbert fc Reeslll988l ). The inter- 
action timescale for these two processes is approximately 
the light crossing time between the two emitting regions, 
and will be in the region of minutes to days for black hole 
masses M BH ~ 10 6 - 1O 8 M . 

Both emission regions are also likely to be correlated at 
some level if they are both modulated by their local accretion 
rate, which varies as accretion rate fluctuat ions propagate 
through the flow (e.g. lArevalo fc Uttlevl2006l ). In a standard 
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accretion disk the timescale for propagation of fluctuations 
between the two emission regions is governed by the vis- 
cous timescale of the disc. For Mbh ~ 10 6 — 1O 8 M this 
times cale will be in the region of weeks to years |Czernvl 
2006). If any of these processes are significant, their effects 
should be apparent from time series analysis of light curves 
from both wavebands. A combination of these processes oc- 
curring at the same time could make individual reprocessing 
models harder to detect. 

Studies of correlations between variations in different 
wavebands are a potentially powerful tool for investigating 
the connecti ons between different emission mechanisms (see 
IUttlevll2006l for a short review). X-ray /optical correlations 
on long timescales have been seen in radio-quiet AGN (e.g . 
Uttlev et all 120031: 1 Arevalo et all I2OO8I ; I Arevalo et al J 120091 ; 



Breedt et alj 20 09). Together with the optical-optical lags 
(e.g. Cackett et al.l 120071 ) they imply that a combination of 



accretion fluctuations and reprocessing produces much of 
the op tical variability. X-ray /UV correlations have been seen 
in e.g. lNandra et all (fl998T ); ICameron et al.l (|2012f ). however, 
there are currently fewer examples than the X-ray /optical 
correlation studies. 

The target of the present paper - the low-mass Nar- 
row line Seyfert 1 (NLS1) galaxy NGC 4051 - has been 
the subject of several such studies. iDone et al.l (|l990l ) ob- 
tained ~ 3 days of contemporaneous X-ray, UV, optical 
and IR data, but found little variability at longer wave- 
lengths despite strong, rapid X-ray variations. Using light 
curves spanning ~ 1 000 days in the optical and X-rays, 



IPeterson et al.1 ( 2000T ) revealed significant optical variabil- 
ity on longer timescales that appeared t o be correlated 
with t he (lo nger timescal e ) X-ra y variations. Sh emmer et al.l 
[|2003h and iBreedt et all (|2010l ) also found a significant X- 
ray /optical correlation, the latter u sing ~ 5000 days of mon- 
itoring data. On shorter timescales iMason et al. I (20021 ) and 
ISmith fc Vaughanl (|2007l ) used ~ 2 day XMM-Newton ob- 
servations to search for X-ray/UV correlations, with incon- 
clusive results. 

The rest of this paper is organised as follows. Section [2] 
discusses the observations of NGC 4051 and extraction of the 
X-ray and UV light curves. Section [3] gives an analysis of the 
variability amplitudes and UV power spectral density, cross- 
correlations are discussed in section [4] and the implications 
of these results are discussed in section [S] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target and observations 

NGC 4051 is a nearby (z = 0.0023) NL S1 galaxy with blac k 
hole mass, M BH « 1.7±0.5 x 1O 6 M jPennev et al .1120091) . 
at a Tully-Fisher distance, D « 15.2 Mpc (|Russel]||2002h . 

Simultaneous, or quasi-simultaneous, observations of 
the two bands are required in order to search for corre- 
lated variability. The XMM-Newton and Swift observatories 
are suited for this task . The thre e EPIC X-ray dete ctors 
fpn; IStriider et al.l l200ll . MPS 1/ 2; iTurner et alj|200ll ) and 
co-aligned Optical Monitor (OM; IMason et aL 2001) makes 
XMM-Newton an ideal instrument for probing the X-ray to 



Table 1. XMM-Nexton observation summary. The columns list 
(1) the spacecraft revolution number, (2) the start date of the 
onservation, (3) the EPIC-pn observation duration, (4) the OM 
observation duration, (5) the number of OM images in each ob- 
servation. 



XMM 


Observation 


EPIC-pn 


OM 


No. OM 


rev. no. 


Date 


On time 


On time 


images 




[Y-M-D] 


W 


N 




1721 


2009-05-03 


45717 


45105 


33 


1722 


2009-05-05 


45645 


45103 


33 


1724 


2009-05-09 


45548 


45003 


33 


1725 


2009-05-11 


45447 


44903 


33 


1727 


2009-05-15 


32644 


31340 


21 


1728 


2009-05-17 


42433 


37367 


26 


1729 


2009-05-19 


41813 


41267 


26 


1730 


2009-05-21 


41936 


40894 


26 


1733 


2009-05-27 


44919 


39168 


22 


1734 


2009-05-29 


43726 


43182 


27 


1736 


2009-06-02 


44946 


44164 


22 


1737 


2009-06-04 


39756 


34574 


23 


1739 


2009-06-08 


43545 


43000 


27 


1740 


2009-06-10 


44453 


43909 


28 


1743 


2009-06-16 


42717 


42116 


26 



UV correlation on short time scale s. The rapid re sponse and 
flexible scheduling of Swift (jGehrels et al.1 12004). with the 



X-ray Telescope (XRT; iBurrows fc The Swift XRT Team! 
20041) and co-aligne d Ultraviolet/Optical Telescope (UVOT; 
Roming et al.ll2005T ). allow for X-ray and UV monitoring on 
longer timescales. 

NGC 4051 was observed by XMM-Newton in 15 sepa- 
rate observations over a period of 45 days during May- June 
2009. Each observation lasted ~ 40 ks giving a total of ~ 580 
ks of usable data. Individual observation details are listed 
in table Q] We make use of the EPIC-pn detector in small 
window mode, with an ~ 4 x 4 arcmin field-of-view. With 
roughly the same start and end time as the OM, this gives 
a total of ~ 530 ks of simultaneous UV and X-ray data. 

All of the OM observations were taken in Imaging Mode 
in the UVW1 filter (central wavelength 2910 A) with 2 x 
2 pixel binning. This gives a field-of-view of 8 x 8 arcmin 
and a spatial resolution of 0.95 arcsec/pixel. The exposure 
length varied between 1000-1400 sec from observation to 
observation, giving ~ 30 exposures per observation and a 
total of 406 across all 15 observations. 

To compliment the XMM-Newton data set, 51 Swift 
ToO observations were made covering the same epoch. The 
observations are separated by ~ 1 day and typically ~ 1.5 
ks long, with 1-3 UVOT exposures per observation, giving a 
total of 71 frames. The UVOT exposures were taken in the 
uvwl filter, which has approximately the same bandpass as 
the OM UVW1, with a field-of-view of 17 x 17 arcmin and 
pixel scale of 0.5 arcsec. The XRT has a 23 x 23 arcmin 
field-of-view with a pixel scale of 0.236 arcsec and covers 
the energy range 0.2-10.0 keV, similar to EPIC-pn. 



2.2 OM light curves 

The Observation Data Files (ODFs) for our target were 
extracted from the XMM-Newton Science Archive (XSA) 
and processed using XMM-Newton Science Analysis System 
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(SAS vll.O.O) routine omichain. Custom made idiQ scripts 
were made to perform source photometry and apply instru- 
mental corrections. Source counts were extracted in a 6 arc- 
sec radius aperture for the galaxy nucleus and 3 field stars 
present in the images. Background counts were extracted in 
a 30 arc second radius aperture, placed in a region away from 
the host galaxy and field stars. Accurate count rates from 
aperture photometry of the OM images can only be pro- 
duced once five ins trumental corrections have been applied 
l|Mason et al.ll200lh . These are, in order of application to the 
extracted counts: the point spread function (PSF1), coinci- 
dence loss (CL), CCD dead-time (DT), the UV point spread 
function (PSF2) and time-dependent sensitivity degrada- 
tion (TDS) corrections. The concatenated background sub- 
tracted light curves for the XMM OM sources and back- 
ground region are shown in Fig. [1] 

In the UVW1 filter there will be a significant contribu- 
tion to the observed nuclear light from the host galaxy. This 
should be constant (to within the random and systematic 
errors of the aperture photometry) and so we have not tried 
to remove it, but as such it should not affect the PSD or 
correlation analysis in any important way. 

As a test of the background subtraction and photome- 
try procedure we tested for (zero lag) correlations between 
the background subtracted light curves of the sources and 
a second background region for all 15 observations. For 
~ 400 data points a Pearson linear correlation coefficient 
r 0.13 indicates a weak but statistically significant corre- 
lation (p < 0.01). We find no significant correlation between 
each of the background subtracted sources, but in the source 
vs background tests, values of r up to 0.3 are observed. The 
strength of this correlation is also observed to change be- 
tween 0.0 and 0.3 when using a different background region. 
The mean correlation coefficient between source vs source 
and source vs background for individual observations is very 
low (—0.1 < r < 0.1). This indicates that the correlation is 
caused by changes in the background over the course of the 
15 observations. We are cautious of this fact during the rest 
of the analysis and use source light curves subtracted us- 
ing various background regions. We find that the choice of 
background region has no effect on any subsequent analysis. 



2.3 EPIC-pn light curves and spectra 

The EPIC-pn ligh t curves used in this a nalysis are the same 
as those used in IVaughan et~al]|201li The raw EPIC-pn 
data were processed from the ODFs using the SAS (vll.O.O). 
Events lists were filtered using pattern 0-4, FLAG = and 
were visually inspected for background flaring. Light curves 
were extracted from the filtered events files using a source 
aperture of radius 35 arcsec and a non-overlapping larger 
background region on the same chip. Light curves were ex- 
tracted with bin size of At = 5 s and an energy range of 0.2- 
10.0 keV. The background subtracted EPIC-pn light curves 
are shown in the bottom panel of Fig. [3] Spectra were ex- 
tracted and binned to a minimum of 25 counts per bin. Re- 
sponse files were created using rmfgen vl.55.2 and afrgen 
v.1.77.4. 
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Figure 1. Concatenated light curves for the XMM OM sources. 
From top to bottom panel they are central nucleus, background 
region and field stars 1, 2 and 3 respectively. 



2.4 UVOT and XRT light curves 

Visual inspection of the 71 UVOT exposures revealed con- 
siderable target movement in 40, leaving 31 usable frames. 
The exposure time varied slightly across the observations 
but is typically ~ 500 s. The HEAsoft (v6.10) tasks uvot- 
SOURCE was used to extract source counts from the galaxy 
nucleus and the same field stars used in the OM photom- 
etry. The optimum extraction ra dius is 12.5 pixel s (~ 6 
arc seconds) for the uvwl filter (jPoole et al . 2008). A 60 
pixel radius background aperture fixed in sky co-ordinates 
in a blank region of the sky was also extracted. Similar to 
the OM, instrumental corrections — PSF1, CL, DT, PSF2 
and TDS — must be applied to the counts extracted in 
each aperture on the UVOT images. UVOTSOURCE automat- 
ically performs background subtraction and applies these 
corrections based on the most up-to-date instrument calibra- 
tion database (CALDB 4.1.2). The background subtracted 
UVOT light curves are shown in the lower panel of Fig. [3] 
We extract XRT counts in th e 0.2-10 keV rane e using the 
online XRT Products Builder (jEvans et al]l2003 ). This per- 
forms all the necessary processing and provides background 
subtracted light curves. 
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Figure 2. Within-observation X-ray and UV variability from the 15 XMM- Newton observations. To better represent the light curves the 
OM (thick black) has been scaled to some arbitrary value and shifted to a mean of 1.5 ct/s. The EPIC-pn (thin red) has been normalised 
to a mean rate of 1 ct/s. The circle and diamond markers represent the mean observation count rate, given by the right ordinate, for the 
OM and EPIC-pn respectively. The four "good" OM observations are indicated by the asteriks next to the revolution number. 
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Figure 3. Between-observation variability for UV (top) and X-ray (bottom) from XMM and Swift. In the top panel the black lines are 
the OM data and the red circles are the UVOT. In the bottom panel the black lines are the EPIC-pn data and the blue dots the XRT. 
In both panels the Swift count rates have been scaled up to account for different effective areas (see section [4.21 for details). 
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3 DATA ANALYSIS 

3.1 Quantifying variability 

In an observed light curve some of the total variance will be 
intrinsic to the source and some will come from variati ons in 
the measurement uncertainties l|Vaughan et al.l [20 031). The 
difference between the two is the 'excess variance', and can 
be used to estimate the intrinsic source variance. The values 
of the variability estimators for the OM and UVOT sources 
are listed in table[2] The variability statistics calculated over 
the whole ~ 45 day XMM-Newton observation show there 
is significant variability in NGC 4051, with F var ~ 8%. On 
timescales within each observation, the variability is weaker 
than on long timescales. 

In 11 of the XMM-Newton observations the excess vari- 
ance in the nucleus and star 1 values are almost identical. 
This indicates that there is a "floor" in the excess variance, 
that isn't accounted for by the errors. We refer to this as a 
sytematic error but are unable to account for it. In the 11 
observations this systematic error is too large to detect any 
intrinsic variability from the nucleus. The variability estima- 
tors in table [5] show the nucleus in the remaining 4 "good" 
observations clearly posses significant variability compared 
to the field stars, and as such will be treated separately form 
the the "poor" 11 observations in the correlation analysis 
in section [4] The "good" UV observations are revolutions 
1725; 1729; 1734; and 1740, and are indicated by an asteriks 
in Fig. El 



3.2 The UV Power Spectrum 

The power spectral density (PSD) describes the amount of 
variability power present in the light curve (mean squared 
amplitude) as a function of temporal frequ ency. The PSD 
of the X-ray light curves is discussed by IVaughan et al.l 
(2011). The UV power spectrum was estimated from the 
15 individual OM observations using standard methods (e.g 
Ivan der Khlll989l ). A 30 ks segment (equal to the shortest 
observation length) was taken from each observation. Within 
each XMM-Newton observation the individual OM expo- 
sures are approximately evenly sampled in time, although 
the exposure times do differ between observations (from 
1200 to 1500 s). The basic periodogram requires evenly sam- 
pled data, and so we interpolated all OM data onto a grid 
evenly sampled at At — 1500 s - the smoothness of the OM 
light curves means that linear interpolation should not af- 
fect the shape of the time series in any significant manner. 
The observed power spectrum may be distorted by leakage 
of power f rom low frequencie s to higher frequencies (van der 
Klis 1989; IUttlev et al.ll2002l ; Vaughan et al. 2003). This can 
bias the data such that the observed spectrum resembles an 
a — 2 power law even if the true power spectrum is some- 
what steeper. 

Figure [4] shows the power spectrum for the NGC 4051 
nucleus, star 1, star 3 and the background for the 15 XMM- 
Newton OM observations. Periodograms were computed 
with absolute normalisation and the Poisson noise level is 
estimated using the formula in Vaughan et al. (2003); Ap- 
pendix A. The background power spectrum is computed for 
the light curve from the background region (the same back- 
ground region used in the source background subtraction) 



Table 2. Quantifying source UV variability. The columns give 
the following information: (1) source name, (2) mean count rate, 
(3) sample standard deviation, (4) excess variance, (5) fractional, 
excess rms amplitude (_F V ar ) . The upper table shows the values ob- 
tained by combining all 15 XMM-Newton OM observations. The 
second and third tables show the values obtained by averaging the 
results from each of the 4 "good" and 11 "poor" XMM-Newton 
OM observations respectively. The lower table shows the results 
from the Swift UVOT observations. 



Object Mean rate a cr^g F va r 

ct/s ct/s [ct/s] 2 percent 



XMM-Newton total 



Nucleus 


40.5 


3.19 


10.1 


7.9 


Star 1 


86.2 


0.55 


0.24 


0.6 


Star 2 


1.6 


0.05 


0.0005 


1.5 


Star 3 


4.6 


0.09 


0.003 


1.2 


Background 


0.4 


0.14 


0.015 


2.6 


XMM-Newton "good" 


4 observation averaj 


res 


Nucleus 


40.0 


0.84 


0.69 


2.1 


Star 1 


86.3 


0.37 


0.09 


0.3 


Star 2 


1.6 


0.05 


0.0005 


0.2 


Star 3 


4.7 


0.08 


0.001 


0.6 


Background 


0.4 


0.006 


0.005 


1.0 


XMM-Newton "poor" 


11 observation avera 


ges 


Nucleus 


40.7 


0.40 


0.14 


0.9 


Star 1 


86.2 


0.42 


0.13 


0.1 


Star 2 


1.6 


0.05 


0.0004 


0.1 


Star 3 


4.6 


0.08 


0.003 


0.9 


Background 


0.4 


0.007 


0.007 


1.8 




Swift total 






Nucleus 


40.1 


4.39 


19.2 


10.8 


Star 1 


48.5 


1.06 


1.01 


0.02 


Star 2 


1.1 


0.04 


-0.001 


0.02 


Star 3 


3.1 


0.09 


0.002 


0.01 



subtracted by a second background region on the opposite 
side of the CCD. In all sources, some power above the noise 
level is present. This is most likely the result of the back- 
ground subtraction issues described in section 12.21 Star 1 
shows a similar red-noise slope, albeit with less power, to 
the nucleus. A cross-correlation test (see section [4| between 
the nucleus and star 1 revealed no significant correlation 
between the two sources. This indicates that the variations 
in star 1 are either intrinsic to the star or caused by the 
problems in the background subtraction. 

Figure[S]shows the resulting NGC 4051 power spectrum 
from the "good" and "poor" OM data, and a two compo- 
nent model is fit to the "good" data. Periodograms were 
computed with fractional rms normalisation (Vaughan et 
al. 2003; Appendix A). We did not subtract the expected 
contribution from Poisson noise but instead included this in 
the model fitting. The simple model comprises a power law 
plus constant to account for the Poisson fluctuations in the 
count rate: P(y) = Av~ a + Pn (where v is the temporal fre- 
quency, A is a normalisation term, a is the power law index 
and Pjv is the power density due to Pois son noise). Th is was 
fitted to the data using xspec v12.6.0 (|Arnaudl fl996). The 
Pn level was allowed to vary freely. Using a y 2 statistic the 
best fit to the data is found to have a = 2.62 ± 0.48 and 
Pjv = 0.059 ± 0.023, with x 2 = 3.7 for 7 degrees of freedom 
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Figure 4. UV power spectrum in absolute units for the back- 
ground subtracted OM sources, produced from the 15 XMM ob- 
servations. The circles, squares, diamonds represent NGC 4051, 
star 1 and star 3 respectively. The triangles are the power spec- 
trum estimate for the background light curves subtracted by a 
background region on the opposite side of the CCD, see text for 
details. The dotted lines are the Poisson noise estimates for each 
source. 
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Figure 5. UV power spectrum produced from the XMM OM ob- 
servations. The black circles are the average of the 4 "good" peri- 
odograms and the grey squares are the average of the 11 "poor" 
periodograms, prior to subtraction of Poisson noise. The red solid 
line is the power law fit to the "good" data. The upper solid black 
line is the Poi sson noise subtracted X-ray (0.2 — lOkeV) power 
spectrum from lVaughan et al.ll201l] . 



(dof). The Poisso n noise level can also b e estimated from the 
formula given in IVaughan et~aH (|2003h which we compute 
to be Pn = 0.05, in agreement with the value derived from 
the PSD. As expected, fitting the model assuming a fixed 
Pjv = 0.05 value gave consistent results for the index param- 
eter (a = 2.50 ±0.35), with x 2 = 3.9 for 8 dof. Errors on the 
model parameters correspond to a 90 per cent confidence 
level for each interesting parameter (i.e. a Ay 2 = 2.7 crite- 
rion). For comparison, the X-ray PSD from Vaughan ct al. 
l|201ll ) is plotted in Fig. [5] 



4 CORRELATION ANALYSIS 

In this section we discuss tests for X-ray/UV correla- 
tions on short timescales within each observation (within 
-observations) and longer timescale (between-observations) . 
Treated individually, each of the 15 XMM observations al- 
lowed us to probe time scales of ~ 2 — 40 ks. Combining the 
XMM and Swift light curves allowed us to look for any long 
term trends in correlation between the two bands over the 
~ 45 days. 

4.1 Within-observation correlations 

Standard time series analysis methods (c.f iBox fc Jenkins! 
ll97d ; |PriestIevlll98ll) require the two light curves to be si- 
multaneous and evenly sampled. This requirement is com- 
plicated by the the OM and EPIC-pn not always starting 
and ending at the same time, the irregular sampling of the 
OM, the read out time of the OM CCD, and any bad OM 
exposures. Where the two light curves are simultaneous we 
linearly interpolate the OM onto an uniformly- sampled reg- 
ular grid. Given that the OM light curves vary smoothly 
within each observation, linear interpolation should not have 
a significant effect on the intrinsic variability. The EPIC-pn 
X-ray counts are re-binned to be contiguous and simulta- 
neous with the OM bins by taking the average count rate 
within the new bin width. We chose a bin width of 1500 s to 
be consistent with the mean sampling rate of the OM (1502 
s) across the 15 XMM observations. The simultaneous light 
curve lengths range from 28.5-43.5 ks. 

4-1.1 The Correlation Function 

The cross-correlation function (CCF) is a standard tool 
for measuring the degree of correlation between two evenly 
sampled time series (x t ,yt) as a f unction of time-lag (c.f 
IBox fc Jenkindll97d ; |Priestlevlll98ll ). We estimated the CCF 
for each XMM-Newton observation individually using the 
idl function C_CORRELATE, shown in panel a of Fig. [6] where 
a positive lag in the plot indicates the UV variations are lag- 
ging those of the X-rays. A large spread in the CCF value for 
any computed time lag is seen. Panel b of Fig. [6] shows the 
average CCF for all 15 XMM observations. The strongest 
feature is the peak in the CCF around ~ 4.5 ks, although 
with a correlation of ~ 0.1, which falls within the confidence 
intervals. Confidence intervals on the average CCF are esti- 
mated using Monte Carlo simulations, where the 95% and 
99% confidence intervals are shown in Fig. [6] These show the 
expected range of CCF values under the assumption that the 
X-ray and UV processes are independent, i.e. in the absence 
of a real correlation. The full details of these simulations is 
given in AppendixfS] The error on the average CCF is given 
by the standard error for N observations at each time-lag t. 

Panel c in Fig. [6] shows the average CCF plot for the 4 
"good" observations (see section [3}- A distinct broad peak 
can be seen around ~ 3 ks with a correlation of ~ 0.5, which 
lies outside the 99% confidence interval. Confidence intervals 
are calculated the same as above except 4 simulated light 
curves are averaged over in each CCF estimate. The small 
error bars on the averaged "good" CCF shows there is little 
scatter in the individual CCF measurements. 

A correlation between optical light curves and X-ray 
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Figure 6. The Cross-correlation function (CCF) for the 15 XMM 
observations. The top panel is the individual CCF for each ob- 
servation, the middle panel is the average CCF of all 15 obser- 
vations, whilst the bottom panel is the average CCF for the 4 
"good" observations. The dashed and dot-dashed lines are the 
95% and 99% confidence intervals respectively, calculated indi- 
vidually for each data subset. Conficence intervals were estimated 
using Monte Carlo simulations of light curves assuming no corre- 
lation (see Appendix lAl 



photon index has been detected in some sources (e.g. 
iNandra et ail [2000), despite there being a weak correlation 
between the optical and X-ray light curves. We therefore 
cross-correlated the UV light curves with the 0.7-2/2-10 
keV hardness ratio (a proxy for photon index) but find a 
CCF shape similar to that between the UV and X-ray light 
curves. This is most likely due to the X-ray spectral shape 
changes being strongly correlated with the overall X-ray 
flux. 



4-1-2 The Cross-Spectrum 

The cross-spectrum i s the Fourier t ransform of the CCF 
(|Box fc Jenkinslll976l , |Priestlevlll98ll '). and has been widely 
used in an alysing X-ray ligh t curves from X-ray bi- 



: ay bi- 
_J Il999t 

iMivamoto fc Kitamotd|l989f ) and more recently from AGN 



naries (e.g Vaughan &; Nowakl Il997l ; iNowak et al 



(e.g. iFabian et al.l 120091 ). It contains the same information 



as the CCF but represents the time-lags and strength of 
correlation in terms of phase difference and coherence as 
a function of temporal (Fourier) frequency. The phase lag 
A(f> can be expressed as a time-lag at a given frequency v: 
r = A0/27O/. Under quite general conditions the phase delay 
estimates are approximately independent at each frequency; 
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Figure 7. The Coherence and Phase-lag parts of the Cross- 
spectrum for the combined 4 "good" XMM observations. The 
dashed and dot-dashed lines are the 95% and 99% confidence 
intervals respectively. Conficence intervals were estimated using 
Monte Carlo simulations of light curves assuming no correlation 
(see Appendix [At 



by contrast, adjacent values of the CCF tend to be corre- 
lated due to the autocorrelation of the individual time series. 
Here we have estimated the cross spectrum using the "good" 
OM data, except that the segments have been trimmed to 
equal length (28.5 ks), corresponding to the shortest simul- 
taneous light curve. The resulting coherence and phase parts 
of the cross-spectrum are shown in Fig. [7] Errors were est i- 
mated using standard formulae (| Vaughan fc Nowakl Il997r i ■ 
and confidence intervals were estimated using simulated 
light curves (see appendix lAl for details). 

The coherence between the two bands is found to be 
low 0.2) at all frequencies. The average time-lag of the 
lowest 5 frequency bins is ~ 3 ks, consistent with what is 
seen in the CCF. A low coherence means that the errors 
on the time delay estimates are most likely underestimated 
using standard formulae, which can increase the apparent 
significance of lags when the errors are estimated using the 
standard formula and the intrinsic coherence is very low 
(e.g. Bendat & Piersol 1986). The cross spectrum is also 
computed for the combined 15 XMM observations which 
gives a coherence consistent with zero for each frequency bin, 
and the average time-lag in the lowest 5 bins is consistent 
with that found with the 4 "good" observations. 



4-1.3 Pre-processing the light curves 

The OM light curves tend to be dominated by slow, 
quasi-linear trends, and these can affect the CCF estima- 
tion (Welsh 1999). We have repeated the CCF and cross- 
spectrum analysis after 'end-matching' the OM light curves 
(i.e. removing a lin ear trend suc h that the first and last 
points are level - see lFougerdli"985h . This 'end-matching' re- 
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moves, to a large extent, linear trends from the data, and 
alleviates the problem caused by circularity of the Fourier 
transform when estimating the cross spectrum. 

We end-matched the XMM OM light curves individu- 
ally and computed the CCF and cross-spectrum using the 
same X-ray light curves as before. The CCF for both the 
"good" and all the data remains mostly unchanged. In the 
cross-spectrum the phase-lag follows the same distribution 
and the coherence remains low in both cases. As this reanal- 
ysis did not substantially alter the results we do not show 
the CCF and cross-spectrum plots here. 

4.2 Between-observation correlations 

Given the extended period and sampling of the XMM and 
Swift observations, we are able to search for possible correla- 
tions and time lags on longer time scales. The XMM-Newton 
data (EPIC-pn and OM) and Swift data (XRT and UVOT) 
were first treated separately, then combined to produce one 
X-ray and one UV light curve. In either case the time sam- 
pling between observations is high ly uneven, and so the Di s- 
crete Correlation Function fDCF: lEdelson fc Krolikl l|l988h ) 
was used to estimate the CCF. We use all the XMM-Newton 
OM data in this part of the analysis, as the variations within 
each OM observation will have little affect on the DCF. 

For the XMM dataset, the midpoint of each original 
OM exposure was used. The EPIC-pn data were binned to 
be contiguous from the start of each revolution, with a bin 
size of 1500 s to be consistent with the mean OM sampling 
rate. A 10 ks DCF bin width is adopted to be consistent with 
the mean OM sampling rate over the extended observation. 
Although the OM exposure length varied between 1200- 
1500 s from revolution to revolution, we treat the source 
count rate in each exposure as a representative of the aver- 
age count rate. As the source varies smoothly in the UVW1 
filter we do not expect this to have any effect on the shape 
of the DCF. We test this by computing the DCF using the 
OM data that was binned onto a 1500 s even grid, and find 
no change in the shape of the DCF. The DCF for the XMM 
dataset for the range —40 < lag < +40 days is shown in 
Fig. [8] where a positive lag means the UV are lagging the 
X-rays. Some peaks can be seen in the DCF but all lie within 
the confidence intervals. The peaks are most likely the re- 
sult of the DCF binning used, combined with the underlying 
shape of the uncorrelated red-noise light curves. 

For the Swift dataset, the UVOT exposures from each 
snapshot were used to represent the mean source count rate 
in the middle of each exposure bin. Again the exposure 
lengths varied from ~ 300-800 s with a mean of ~ 500 s, 
but due to the steepness of the red-noise power spectrum 
this will have no effect on the shape of the DCF as long 
as the DCF bin size is much greater than the mean UVOT 
exposure length. The XRT counts are taken from each snap- 
shot and have a mean exposure length of ~ 500 s. The DCF 
for the Swift data is plotted in Fig. [8] The plotted 95% and 
99% confidence intervals are calculated using simulated light 
curves following the method outlined in appendix [Al 

To make the most of the observational coverage we com- 
bine the XMM and Swift datasets and recompute the DCF. 
As the effective areas of the UV and X-ray instruments on 
either telescope are not identical the count rates from one 
telescope need to be scaled before the DCF can be com- 
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Figure 8. The discrete correlation function (DCF) for the 
between-observations (solid black line). The dotted black line is 
the DCF for the end-matched data. The dashed and dot-dashed 
lines represent the 90% and 95% confidence intervals respectively. 



puted. We estimate this scale factor using the 3 occasions 
the observations overlap and find OM ~ l.lxUVOT, and 
EPIC-pn 15 x XRT. The scaling factor for the X-ray cam- 
eras is consistent with that calculated by WebPIMMlfl No 
reference for a scaling factor between the UV cameras could 
be found, but we find the choice of scaling factor within the 
range ~ 0.5 — 1.5 has no effect on the shape of the DCF. 

In Fig. [3] it can be seen that the UV light curves show 
a gradual increase in counts over the extended observation. 
To account for any underlying long-term trends in the UV 
variability we 'end-match' the overall UV light curve and re- 
compute the DCF for the individual and combined datasets. 
These are shown as the dotted black lines in Fig. [5] 



5 DISCUSSION 

Using UV and X-ray data from XMM-Newton and Swift we 
have analysed the light curves of NLS1 galaxy NGC 4051 
to search for correlations in the variability between the two 
bands. UV variability is detected on short and long time 
scales, however the fractional rms amplitude is smaller than 
that in the X-rays. On days-weeks timescales the fractional 
variability of the UV is F vai ~ 8%, and on short (~ hours) 
timescales F^^ ~ 2% (from the "good" OM observations). 

The excess variance in 4 of the 15 XMM-Newton OM 
observations is found to be considerably greater than the 
remaining 11. The "poor" 11 observations show there is a 
"floor" to the excess variance that isn't accounted for in the 
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errors. In the "poor" observations any intrinsic source vari- 
ations are masked by the errors, and inclusion of these ob- 
servations will weaken the detection of any correlated emis- 
sion. Although 4 out of the 15 observations is a relatively 
small subset, we find the variability statistics of the 4 "good" 
observations clearly very different from the remaining 11 
"poor". The similarity in the overall shape of the CCF for 
the 4 "good" observations is hard to explain as arising by 
chance if the they were all representative of uncorrelated 
emission. Nevertheless, the interpretation of the time lag is 
still treated with some caution. 

Analysis of the UV power spectral density reveals a 
red-noise light curve with a power-law slope of index a — 
2.62 ± 0.48 for the 4 "good" OM observations. We searched 
for correlations between the two bands on time scales up to 
~ 40 ks, treating all the XMM-Newton and the 4 "good" 
observations separately. The CCF for the 4 good observa- 
tions revealed a significant peak of ~ 0.5 at a lag of ~ 3 
ks. Using all 15 XMM-Newton observations the CCF re- 
vealed a weak correlation (~ 0.1) with a peak at ~ 3 ks. 
The cross-spectrum showed the lowest 5 frequency bins to 
have a low mean coherence of ~ 0.2 and a mean phase-lag of 
~ 3 ks in both cases. Combining the XMM-Newton and Swift 
datasets we searched for correlated emission on timescales 
up to 40 days and find no significant correlations. The lag 
significance in the above results was estimated using simu- 
lated light curves. A correlation coefficient of r ~ 0.5 means 
the amount of UV variance prdictable from the X-ray var- 
ince is r 2 ~ 0.25. As the coherence is a "square" quantity, 
this value is consistent with the ~ 0.2 form the coherence. 

From a ~ 12 yr monitoring cam paign on NGC 4051 
using ground based optical photometry iBreedt et all (|2010t ) 
estimated the PSD in the frequency range ~ 10" 8 -10~ 3 Hz. 
In their Fig. 4 they fit an unbro ken power law to th e PSD 
with a = 1.4±S;a. In the paper IBreedt et al.l (|201Ch fitted 
a single-bend p ower-law model to the PSD, as is observe d 
for X-rays (e.g. iMcHardv et"al]|2004 IVaughan et afll201ll ). 
Whilst they do not rule out the single-bend model, they 
find their data is more consistent with an unbroken power 
law. Their Fig 5 shows the acceptance probabilities for the 
single-bend model as a function of high frequency slope an 
and bend frequency vb- Taking our value of a ~ 2.5 as 
the high frequency slope this would give a break frequency 
vb ~ 10 -6 ' 5 Hz. Our PSD is better constrained in the high- 
frequency range (~ 10~ 5 — 10 _3 Hz) and so the slope value 
ct — 2.62 ± 0.48 is consistent with their single-bend model. 

In a samp le of 4 AGN using Kepler data, 

iMushotzkv et ail 1201 ll estimated power spectral slopes 
(assuming a single power law model) of ~ 2.6-3.3 to the 
optical PSD in the ~ 10" 6 5 - 10" 3 5 Hz range. They do 
not attempt to fit a single-bend power-law model to their 
PSDs, but the break frequency for the larger black hole 
masses (~ 10 7 Mq ) in their sample would likely occur at 
lower frequencies than they estimate in the PSD. 

The X-ray PSD in Fig. [5] shows orders of magnitude 
more variability power than the UV. This is consistent with 
what is seen in optical PSDs, where the high-frequency 
power is much less than in the X-rays, although the am- 
plitudes can be s imilar (or even grea ter) at low frequencies 
(e.g, NGC 3783, lArevalo et al]|2009h . A study with simul- 
taneous UV and X-ray coverage on longer timescales is still 
lacking. A bend can be seen in the X-ray PSD at ~ 2 x 10~ 4 
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Table 3. UV and X-ray rms luminosity for the 15 XMM-Newton 
observations. 



Band 


L 

10 41 erg/s 


^rms 

10 41 erg/s 


UVW1 


3.4 


0.3 


0.2-10 keV X-ray 


7.3 


3.9 


1-1.2 keV X-ray 


0.3 


0.2 


5-6 keV X-ray 


0.5 


0.2 



Hz <|Vaughan et al.ll201lf ). If a break was present in the UV 
PSD, it would be seen to occur at much lower frequencies 
that of the X-ray due to the radius of UV emission being 
much greater than that of the X-rays. 

To assess whether the observed X-ray variations are sig- 
nificant enough to produce the variations seen in the UV 
band, we compare the root-mean-square luminosity varia- 
tions in both bands. If the luminosity variations in the UV 
band are greater than the luminosity variations in the 0.2- 
10.0 keV band, then this would in effect rule out the 0.2-10.0 
keV X-ray variations being the dominant cause of variations 
in the UV band. The values in table [3] show the integrated 
X-ray luminosity is greater than in the UVW1 band, and 
the X-ray luminosity variations are a factor ~ 10 greater. 
It is worth noting here that the UVW1 filter is very narrow 
compared to the X-rays. The X-ray band covers a factor of 
~ 50 in wavelength, the UVW1 band covers only a factor 1.3. 
This largely explains the apparently low luminosity in the 
UWV1 compared to X-rays. The ratio of the FWHM to the 
central wavelength of the UVW1 filter is 620A/2910A« 0.2. 
Table [3] gives the rms luminosity for X-ray bands of compa- 
rable fractional energy range to the UVW1 filter. The rms 
luminosity in the narrower X-ray bands is now comparable 
to that of the UV band, albeit with lower mean luminos- 
ity. This shows that, in principle, the X-ray variations could 
drive variations in the UV band. 

Given th e published black hole mass (A/bh ~ 1.7 ± 
0.5 x 10 6 Mq iDennev et alj|2009h it is possible to make pre- 
dicted lag estimates for each reprocessing scenarios based 
on standard disc equations to find the distance of the UV 
emitting. In the Compton up-scattering scenario the lags 
can be expected to be seen in the ~ 1.5-7 ks range for as- 
sumed accretion rate as a fraction of Eddington of 0.01-0.1. 
In the thermal reprocessing scenario the time-lags depend 
on the luminosity of the X-ray band and are expected to 
be ~ 7 ks. The direction and magnitude of our lag from 
the "good" data is consistent with the thermal reprocessing 
scenario. Although the expected time delay of ~ 7 ks is pre- 
dicted from the toy model, the model assumes that the disc 
is heated solely from the incident X-rays, which are them- 
selves coming from a radius r = 0. Both these assumptions 
are not likely to be true for a real AGN. An extended corona 
will increase Rx and a viscously heated disc will decrease 
Ruv and hence the light travel time between the two emit- 
ting re gions. In the propagati ng accretion rate fluctuation 
model l|Arevalo fc Uttley 2006) the timescale of mass flow is 
dictated by the viscous timescale. This is dependent on the 
assumed viscosity parameter and scale hei ght of the geo met- 
rically thin, optically thick accretion disc (|Czerny||2006h . We 
estimate this to be in the region of ~ weeks — years. 

Given the quality of the UV data in the 4 "good" ob- 
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servations, a lag in the region of ~ 1.5-7 ks would have 
manifested itself in the cross-correlation analysis. If the lag 
estimate from the 4 "good" observations is to be believed, 
then crudely speaking ~ 25% of the UV and X-ray variance 
are correlated on times cales of days. This is consistent with 
the lBreedt et all (|201Ch result, where an optical — X-ray cor- 
relation of ~ 30% is reported on timescales of ~ weeks. 
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then selected at random from this light curve and points 
then sampled from this coinciding with times and lengths of 
the real XMM and Swift observations. 

In order to set an upper limit on coherence we simu- 
lated light curves with added variance i.e. some fraction A 
of the simulated X-rays was added to the simulated UV. 
The variances were normalised before the two light curves 
were added. We simulated 10 4 light curves in each band for 
a range added variance and recorded the mean coherence 
value for the lowest 5 Fourier fr equencies. The bias in co- 
herence |Bendat fc Piersolll986l . section 9.2.3) is given by 
B[j 2 ] = l/n(l- 7 2 ) 2 where 7 2 is the coherence and n is the 
number of segments going into the cross-spectrum (15 in our 
simulations). When the coherence is low the bias dominates 
and acts to shift up the observed coherence value and must 
be subtracted from the computed value. The distribution in 
coherence values for each A were then compared to our ob- 
served coherence value ■Jobs- The mean coherence value for 
the distribution of coherences where 90% of the values fall 
above 7 2 6s is taken as the upper limit. 



APPENDIX A: SIMULATING LIGHT CURVES 

The red-noise nature of UV light curves means that individ- 
ual data points in the light curve are correlated with adja- 
cent points. Confidence intervals were placed on the CCF 
and cross-spectrum measurements from the 15 XMM obser- 
vations using Monte Carlo simu lations of uncorrec ted light 
curves. We used the method of iTimnier fc Konid (|!995l 1 to 
simulate 10 4 light curves in each band with length 50 days, 
the same time resolution as the re-binned data (At = 1500s) 
and appropriate PSD shapes. The X-ray PSD was modelled 
by a bending power-law with low frequency slope -1.1, high 
frequency slo pe -2.0 and break freq uency 2 x 10~ 4 Hz pa- 
rameters from Vaughan et al.l i 20 111 ) . The observed rms-flux 



relation (see lUttlev et all 2005 ) was added to the simulated 



X-ray light cu rves by computing the exponential functio n 
of each point (|Uttlev et al.ll2005l ; IVaughan fc Uttlevl F2008) 



The UV PSD was modelled with an unbroken power law 
with slope -2.1 (see section 3). Observational noise was 
added to each simulated UV and X-ray light curve by draw- 
ing a Poisson random deviate with mean equal to the mean 
count per bin in the real light curves. We then took 15 seg- 
ments corresponding to the times and length of each real ob- 
servation from the 50 day simulated light curves. The CCF 
was computed for each segment before averaging. This re- 
sults in 10 4 simulations of the averaged CCFs from which 
we extracted the confidence intervals at each lag. Confidence 
intervals on the cross-spectrum lag and coherence were com- 
puted using the same approach. 

Light curves were simulated to put confidence intervals 
on the DCF using the above procedure, except the generated 
light curves were 20 times longer than before (1000 days). 
The time resolution was 1500 s for XMM and 500 s for Swift 
and 500 s for the combined datasets. A 50 day segment was 



